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0> ' Abstract 

The neutron charge form factor G^{Q'^) is investigated within a constituent quark 
model formulated on the light-front. It is shown that, if the quark initial motion 
^ ' is neglected in the Melosh rotations, the Dirac neutron form factor F^{Q'^) receives 

PJ . a relativistic correction which cancels exactly against the Foldy term in G^{Q'^), as 

it has been recently argued by Isgur. Moreover, at the same level of approximation 
the ratio of the proton to neutron magnetic form factors G^(Q^)/G^((5^) is still 
given by the naive S'[/(6)-symmetry expectation, —3/2. However, it is also shown 
P^ \ that the full Melosh rotations break SU{6) symmetry, giving rise to G^{Q^) ^ and 

G^{Q'^)/G^{Q'^) 7^ —3/2 even when a 5C/(6)-symmetric canonical wave function is 
assumed. It turns out that relativistic effects alone cannot explain simultaneously the 
experimental data on G|(g2) and Gl^{Q^)/Glj{Q^). 



"To appear in Physics Letters B. 



The elastic nucleoli form factors contain important pieces of information on the inter- 
nal structure of the nucleon and therfore an extensive program aimed at their experimental 
investigation is currently undergoing and planned at several existing facilities around the 
world [|l|. In what follows we will focus on the so-called nucleon Sachs form factors, which 
are defined as 

GUQ') = F,\Q') + F,\Q^) (1) 

where FflQ"^) [F^{Q'^)] is the Dirac [Pauli] nucleon form factor, appearing in the usual 
covariant decomposition of the nucleon electromagnetic current matrix elements, viz. 

<mp',s%LmN{p,s)>=u{p',s')[F,'^{Q')r + Fi'{Q')'-^]u{p,s) (2) 

with Q^ = —q ■ q and M being the squared four-momentum transfer and the nucleon mass, 
respectively. As it is well known, the nucleon Sachs form factors may be interpreted in 
the Breit frame as the Fourier transforms of the nucleon charge and magnetisation density, 
respectively. In this respect, from Eq. (|l]) the squared nucleon charge radius, r%, is therefore 
given by the sum of two terms, namely 



r%^-6 



dG^jQ 



Q2=0 ^^^-' 



where rfjy = —6[dFf{Q'^)/dQ'^]Q2=Q and k^ = F^{0) is the nucleon anomalous magnetic 
moment. The second term in the r.h.s of Eq. (|^) is usually referred to as the Foldy contri- 
bution and it is of relativistic origin. 

In case of the neutron the charge radius has been nicely determined from electron- 
neutron elastic scattering experiments at very low energy, obtaining r^ = — 0.113±0.005 fm^ 
0). Since from Eq. @ the Foldy contribution turns out to be ~ —0.126 fm"^, the exper- 
imental value of the neutron charge radius appears to be almost totally explained by its 
Foldy term alone, i.e. by relativistic effects. This result, which implies a small value for 
rfpf (~ 0.013 ± 0.005 fm^), has been viewed [^ as an indication of the smallness of the 
intrinsic charge radius related to the neutron rest-frame charge distribution. Nevertheless, 
very recently [§] the interpretation of the neutron charge radius as arising from its internal 
charge distribution has been asserted again. Indeed, it has been argued that, going beyond 
the non-relativistic limit when the Foldy term firstly appears, the Dirac neutron form fac- 
tor F^{Q'^) receives a relativistic correction that cancels exactly against the Foldy term in 
G'%{Q'^). Such a statement has been inferred from the observation that the well-known phe- 
nomenon of zitterbewegung, which produces the Foldy term, cannot contribute to the charge 
radius of the neutron, because the latter has zero total charge 0. 

The aim of this letter is to address the issue of the relativistic effects on G%{Q'^) 
adopting a constituent quark {CQ) model formulated on the light-front. 



Let us briefly recall the basic notation and the relevant structure of the nucleon wave 
function in the light-front formalism. Following, e.g., Refs. [0, |^, §| the light-front nucleon 
wave function is eigenstate of the non-interacting angular momentum operators j^ and jn, 
where the vector n = (0, 0, 1) defines the spin quantization axis. For a system of three quarks 
with equal masses the squared free-mass operator is given by Mq = J2^=i{kfj_ + m^)/^i, where 
m is the CQ mass, C,i = pt / P^ ^"^^ ^i± = Pi± ~ ^iP± are the intrinsic light-front variables. 
The subscript ± indicates the projection perpendicular to the spin quantization axis and the 
plus component of a 4-vector p = {p^,p) is given hj p'^ = p^ + n ■ p; finally P = {P~^, P±) = 
Pi + P2 + P3 is the light-front nucleon momentum and pi the quark one. In terms of the 
longitudinal momentum kin, related to the variable ^j by kin = [CjMd — {kfj_ + m^)/^jMo] /2, 
the free mass operator acquires a familiar form, viz. Mq = X^^i Jm'^ + kf = J2i=i Ei with 

ki = {ki±,kin). Disregarding the colour degrees of freedom, the light-front nucleon wave 
function can be written as 



({6^^±;^:r,}|vi/^-) 



E1E2E3 



^0^166 {,^^ 



E {{'^lm^\{^^}){{k^■,^^r,}\x7) 



(4) 



where the curly braces { } mean a list of indices corresponding to z = 1,2,3; z/j (ri) is 
the third component of the quark spin (isospin); TV = 11^=1 Rjikj±, Q, rn) is the product of 
individual (generalised) Melosh rotations, viz. 

m + ^jMo — ia^^'i ■ {h x kj±) 



Rjikj±,^j,m) 



m + ^jMoy + B 



(5) 



i-L 



with a being the ordinary Pauli spin matrices. In what follows we will limit ourselves to the 
case of a pure SU{6) symmetric canonical (or equal-time) wave function, namely 



{{k^;l^^Ti}\x''N) = Ws{k,p) 



V2 



$ 



00 



({z/.r,}) + $1;^, r{z/,r,}) 



(6) 



where k = {ki ~ k2)/2 and p = k^ are the Jacobian internal co-ordinates for the three- 
quark system, and ws{k,p) is a completely symmetric radial S-wave function. Finally, the 
spin-isospin function ^^^r^^ {{i^iTi}) is defined as 



Ms ^ ^ z ^ 

A'lj' 



(7) 



where 5*12 {T12) is the total spin (isospin) of the quark pair (1, 2). The normalisation of the 

-f -, -f -f -f -f -f 2 -> -► 

wave function (H) is: J2{,yiTi} I dki dk2dks6{ki+k2+k3) ({^j, z^jTjJIxaT) = J dk dp ws{k,p) 

= 1. 



As in Refs. [0, ^ we consider the one-body component of the electromagnetic (e.m.) 
current operator including CQ form factors, namely 

3 3 r T^^c 






a" 



Lriq 



Qfi 



2m 



fHQ' 



where a'^^ = |[7'^,7^], Gj is the charge of the j-th quark, Kj the corresponding anomalous 

magnetic moment and ff(2){Q'^) its Dirac (Pauli) form factor (with fl(2)i^) ~ ^)- ^^ i^ ^^^^ 
known that in the light-front formalism the form factors for a conserved current can be 
obtained using only the matrix elements of the plus component of the current operator (see, 
e.g., Ref. 0) and, moreover, for Q^ > the choice of a frame where q~^ = allows to 
suppress the contribution of the Z-graph (i.e., pair creation from the vacuum) 0. More 
precisely, in what follows the four-momentum transfer q is taken to be purely transverse, i.e. 
q = (0, q±, 0) with Q'^ = — g^ = |glp. In case of the nucleon one has 



(vi>;-|x+ |vi>5;-)=Fi^(g2)5^ 



V NVN 



^^^2^(<5^) {^'nWvI^n) 



(9) 



and therefore FflQ"^) is related to the non spin-flip matrix elements, while F-^IQ"^) can be 
extracted from the spin-flip ones. Using Eqs. (H-H) the nucleon Dirac and Pauli form factors 



are explicitly given by ||T0| 
Ff("^(Q') 



iJm [dfcj [dk'^] \/|j§§^ y^sik',P') wsiKp) I IIA/'W. 



eu{d). 



~. Q 



Im 



A1A2 + Bi- B2 

3AiA2-Bi-B2] A3 + 2{AiB2 + A2Bi) -Bs 
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+ 



Q 



+ 9 



Fl^^^Q') 



2 [A1B2 + A2B1) ■ B 
3M 



2m 



{3A1A2 - Bi ■ B2) A3 



(10) 
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[d^] [dkA_] [dk'j 



I EiE2E3M'o 



w:,{k',p')ws{k,p} {U-^'jK 



A1A2 + Bi- B2 



fU(d), 



Q 
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2 {B^yB2 + B2yBA ■ B3 - B3y {A1A2 + Bi-B2] + 
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2A3 {A,B2y + A2B,y)] [e„(,)/r^'HQ') + 2erfH/f "Hg2 

2 \BiyB2 + B2yBij ■ B3 — B^y {A1A2 + Bi- B2 
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2A3 {A,B2y + A2B 



is/> 



2m, 



K,, 



fU{d) 



d(u). 



u{d)f2 {Q ) + 2fi;d(u)/2 {Q 



(11) 



where 



m 

dk\ 



dk'i 



dki± dk2± dk3± 6[ki± + k2_L + hi] 

dk'i_i_ dk'2i_ dk'2.1. 6[k'i±_ - ki_L + ^igl] 6[k'2i_ - k2± + 6'f±] 

S[k';^ - h± + (^3 - 1) q±] 



(12) 



and 



AT, = l/^{m + i,M,y + kl 
A, 



jy 



A3 
B 

Bjy 

Bzz 



l/^(m + OM'o)2 + A;'^'^ 
(m + ijM'o) (m + ijMo) + k'j^ ■ kj^ 
(m + ^3M'o)k3x - (m + 6Mo)/c'3x 
(m + ^jMo)k'jy - (m + ^jM'o)k^y 
(m + ^,M'o)k,x - (m + ^jMo)k',x 

k' ■ k- — k' ■ k- 

3^ jy jy j^ 

k 3xk3y + k syksx 

-(m + e3M'o)(m + ^sMo) - k'^^ks. + k'syk^y 

(m + ^3M'o)k3y + (m + ^3Mo)k'3y 



(13) 



In the non-relativistic hmit, assuming both point-hke constituents and the SU(6)- 
symmetric wave function (^, the Fourier transform of the non-relativistic charge density, 
E.Lie, 6if-r-), yields G^iQ') ^ EU^MQ'), i-e. G^(Q') ^ F^iQ') and G^(Q2) -^ 0, 
with 



FciQ' = |gt 



rf/c dpws{k,p + -q) ws{k,p) 



(14) 



In a similar way, in case of the non-relativistic magnetisation density the nucleon magnetic 
form factors are given by Gm(Q^) -^ 3Fc(Q2) and Gl^iQ^) -^ -2Fc{Q^). Note that the 
SU{6) symmetry predicts both G%{Q^) = and Glj{Q^)/Glj{Q^) = -3/2. 

To reach the non-relativistic limit from Eqs. (|10|-|11|) is not a trivial task; moreover, 
as argued in Ref. 0, a subtle cancellation of relativistic corrections is expected to hold 
when the Foldy term firstly appears. We now want to find out the assumptions that allow 
to obtain the specific cancellation claimed in Ref. |0 and, then, we want to carry out the 
appropriate non-relativistic reduction of Eqs. (|Ty-|TT|). The basic assumption made in Ref. 
is to neglect the transverse motion of the constituents in order to avoid spin-fiip effects 
from the Wigner rotations of the constituent spins. Thus, our first assumption is to put 
kj± = in the Melosh rotations (|), yielding for j = 1, 2, 3 



Rj{kj^,Q,m) 



(15) 



As for the Melosh rotations in the final state (i.e., after virtual photon absorption), we 
cannot make the same assumption ([T5|), otherwise in case of point-like constituents only a 
vanishing F-^^Q"^) could be obtained[[ Thus, inspired by Ref. 0, we consider the following 
approximation for j = 1, 2, 3 



R,{k',^,^„m) -. . / „ '^' (16) 



2m — ia'^^^ ■ (fi x fc'j^) 

'i-L 

where the denominator is included to maintain the correct normalisation of the final light- 
front wave function. Note that the term proportional to a^^^ in the r.h.s. of Eq. (|TB|) is 
clearly reminiscent of the non-relativistic magnetisation current. Finally, we have to specify 
the final transverse momenta k' jx appearing in Eq. (|TB|). To this end we neglect the initial 
longitudinal motion, i.e. we consider kjn = 0, which together with the assumption kj± = 
yields C,j = 1/3. Thus, from the delta functions in Eq. (|12D the final transverse momenta turn 
out to be: k'i± = k'2± = — ^±/3 and fc'sj, = 2q±/3. Our assumptions ( p!^ - |T6| ) correspond 
to consider ^jMq = ^jM'q = m in Eq. (||), leading to: Afj = l/2m, Af'i = M'2 = 



l/y'4m2 + QV9, AT'g = l/^4m2 + 4QV9, Aj = Am", A3 = 2mQ/3, B^y = B2y = 2mQ/3, 

B^y = —AmQ/S, B^y = —Am'^, while all the other i?'s and -B's are identically vanishing. 
Finally, for compatibility with the non-relativistic reductions 6*^^(0) —>■ 3 and G2^(0) -^ —2 
we neglect any binding effect in the nucleon mass, i.e. we consider M = 3m (cf. also Refs. 

ii). 

Before carrying out the effects of the assumptions (|T5|-0) on Eqs. ( p!OHTT|) , we have to 
restrict ourselves to the case of S'[/(2)-symmetric CQ form factors, otherwise any deviation 
of G^iQ"^) from zero could be attributed to a possible flavour-dependence of the CQ internal 
composite structure. Thus, in what follows we will consider 

fm') = fiQ') 

kjHQ') = e,Kf{Q') (17) 

where f{Q'^), fiQ"^) and k do not depend on the flavour of the CQ (and /(O) = /(O) = 1). 
Using all the above- discussed assumptions it is straightforward to obtain from Eqs. (|TO|-|TT]) 
the following explicit expressions 

Fim = Fo{Q') + 3K^^FoiQ') 

F-{Q-) = -2F,{Q^)-^^-—-2nF,{Q^-^' 



AM'^ + Q^ "'^ MM2 

F^iQ') = 2Fo{Q')+3kFo{Q') 

4A/f2 

F,-m = -2F,m ^j^2 ^ Q2 - ^'^^om 



A vanishing F2 {Q ) is appropriate in the case of the heavy-quark Umit m ^ 00 (cf. Ref. ]l 1|) . 



where 



FoiQ') = 



f{Q' 



1 + QyM^ 



[d^] [dkj_] [dk'j 



E'^E'^E'^Mo 



w*s{k',p')ws{k,p) (19) 



while Fo(Q^) is given by Eq. (|T^) but with f{Q^) replaced by f{Q^). In terms of the Sachs 
form factors (|1|) one gets 



GW') - 


- (1 


- 


^,„)FM-') 


Gim - 


= 


GlsiQ') - 


= 3 


Fo(Q2) +«:(! + 


Gi,m - 


= — 


2 


Fa(Q'') + ri[l 



4M2 



)h{Q^ 



4M2 



)^o(Q' 



(20) 



The relativistic corrections to Ff2{Q'^), contained in Eq. (|T8[), largely differ for proton 
and neutron, but: i) F^iQ"^) receives a relativistic correction that cancels exactly against 
the Foldy term, -F^{Q'^) ■ Q'^/AM'^, so that the S'f/(6)-symmetry prediction Ge(Q^) = 
still holds, as argued in Ref . ; the same type of cancellation does not occur in case of the 
proton charge form factor G^^iQ"^)] ii) the ratio G^j^Q"^) / G^iQ"^) is still given by the simple 
S'f/(6)-symmetry expectation, —3/2. Both results are independent of any particular choice 
of the CQ form factors (both Dirac and Pauli ones), provided the latter are taken to be 
SU{2) symmetric (see Eq. (|17D). Note that Eq. (|20|) predicts a vanishing value for G^^iQ"^) 
at Q^ = 2M2; this is completely at variance with experimental data and signals that the 



applicability of Eqs. (0) and ( pO]) is limited only to low values of Q^ (i.e., Q^ « M^). 
Finally, note also that in Ref. [^] it is claimed that at first order in Q"^ /M'^ the leading term 
in G%{Q'^) is given by its Foldy term. Such a statement is incorrect, because it is obtained 
using directly in G^{Q'^) the non- relativistic limits of F" and F2{Q'^)i ignoring in this way 
the Isgur's cancellation mechanism of the Foldy term. 

We can now carry out the non-relativistic limit of the nucleon form factors by con- 
sidering the formal limit M — i> oo in Eqs. (plsj) and (^0|). Taking also into account that the 
non-relativistic reduction of the integral in Eq. (p!9| ) leads to Eq. (jl^) (see Ref. |^), one has 



FliQ') - 


- FciQ')fiQ') 


F^Q") - 


-^ 


F',m - 


-. F,iQ')[2fiQ') + 3KfiQ')] 


F^Q') - 


-. ~2F,m [f{Q') + ^fm] 



(21) 



and 



GUQ' 



Fcm f{Q' 



GliQ') - 

G^(Q2) ^ -2F,(Q2) [/(Q2) + /./(Q2)] (22) 

It is now worthwhile to make the following two observations: 



the approximation leading to Eqs. (pTSD and (^) is mainly based on neglecting the 
constituent initial transverse motion in the Melosh rotations. Thus, Eq. (^8|) can 
represent a good approximation of the full calculations (p!0| - pJ] ) only when the average 



value of the transverse momenta, (p±), is much smaller than the constituent mass 
m. However, in QCD both m and {p±) are expected to be of the order of the QCD 
scale, Aqcd ~ 300 MeV . Moreover, in quark potential models {p±) turns out to be 
significantly larger than m, because of the high momentum components generated in 
the light-baryon wave functions by the short-range part of the effective quark-quark 
interaction 0; 

• the Melosh rotations break in general the SU{6) symmetry. Indeed, these rotations, 
being momentum and spin dependent, produce a re-coupling of the constituent orbital 
angular momentum and spin; in other words, even if the canonical wave function 
(^ is assumed to be factorized into a spatial part times a spin-isospin wave function 
(which can be classified according to the SU{6) multiplets), after the application of the 
Melosh rotations (^) the light-front wave function (§) cannot be any more expressed 
as a product of a spatial part times a spin-isospin function. As a result, the light-front 
wave function (^) is not SU{6) symmetric and therefore we do not expect in general 
to have G%iQ') = and G^(Q2)/G^(Q2) = _3/2. 

Consequently, the inclusion of the effects of the quark initial transverse motion in Eqs. ( [T0| - 
0) could lead to ^^(Q^) ^ and Glj{Q'^)/Glj{Q^) ^ -3/2. The former point is clearly 



illustrated in Fig. 1, where the results of the calculations of Eqs. (|lO| - pn] ), performed for 



point-like CQ's adopting a simple gaussian-like ansatz for the radial function ws{k,p) oc 



exp[— (/c^ + 3p'^/4:)/2a'jjQ], are reported for various values of the quantity {p±) = j4:/3 ajjo 
at m = 220 MeV (chosen as in Ref. [0). It can clearly be seen that relativistic effects 
may contribute significantly to G'^{Q'^) when (p^) > -mQ Note that, as {p±) increases, the 
neutron charge radius appears not to exceed ~ 40% of its experimental value. 

The SU{6) breaking associated to the Melosh rotations heavily affects also the ratio 
G1i{Q'^)/Gm{Q^)j ^^ it is illustrated in Fig. 2. It should be reminded that the experi- 
mental data on G^(Q^) and G^{Q^) exhibit the well-known dipole behaviour, leading to 
Glj{Q'^)/Gl^{Q^) ~ Up/ fin ^ -1.46 with only a 10^15% uncertainty up to Q^ _ 1 {GeV/cf 

^This result does not depend on the particular choice of a gaussian-like ansatz for the radial function 
wsikiP), provided the average value {p±) is kept the same. In particular, the nucleon eigenfunction corre- 



sponding to the quark potential model of Ref. |12 , yields {p±) ~ 0.58 GeV; we have checked that the use 
of the S'C/(6)-symmetric part of this wave function leads to results almost coinciding with the solid curve of 
Fig. 1, which corresponds indeed to {p±) ~ 0.58 GeV. 



(cf., e.g., Ref. |I|). Thus, Figs. 1 and 2 clearly indicate that relativistic effects alone cannot 
explain simultaneously the experimental data on G^(Q^) and G^^.j{Q'^)/G^{Q'^); therefore, 
other mechanisms, like, e.g., the presence of the mixed- symmetry S-wave generated by the 
spin-spin forces among CQ's and/or the effects of flavour-dependent CQ form factors and/or 
the contribution of many-body e.m. currents, have to be invoked. 

In conclusion, the neutron charge form factor G%{Q'^) has been investigated within 
a constituent quark model formulated on the light-front. It has been shown that, if the 
quark initial motion is neglected in the Melosh rotations, the Dirac neutron form factor 
F"(Q^) receives a relativistic correction which cancels exactly against the Foldy term in 
G'^(Q^), as it has been recently argued in Ref. p. The same type of cancellation does 
not occur in case of the proton charge form factor G^((5^). Moreover, at the same level of 
approximation the ratio of the proton to neutron magnetic form factors G^(Q^)/G^(Q^) is 
still given by the naive S'?7(6)-symmetry expectation, —3/2. These results are independent 
of the electromagnetic structure of the constituents, provided their form factors are SU{2) 
symmetric. However, since the full Melosh rotations break SU{6) symmetry, both G^IQ"^) ^ 
and G\.^{Q'^) / G^{Q'^) ^ —3/2 can be obtained even when a S't/(6)-symmetric canonical 
wave function is assumed. It has been shown that relativistic effects alone cannot explain 
simultaneously the experimental data on G'^(Q^) and G^.jiQ"^) / G^^iQ"^) . 
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Figure 1. The neutron charge form factor G^iQ"^) versus Q^. The experimental data correspond 
to the results of the analysis of Ref. |13| performed in terms of the Reid-Soft-Core nucleon-nucleon 



interaction. The various curves are the results of the calculations of Eqs. (10-11), obtained assuming 
point-like CQ's and adopting a simple gaussian-like ansatz for the radial function ws{k,p) (see text). 
The value of the constituent quark mass has been chosen to be m = 220 MeV from Ref. [l^. The 
dot-dashed, dotted, dashed and solid lines correspond to {p±) = ^4/3 ano = 0.23,0.35,0.46 and 
0.58 GeV, respectively. 
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Figure 2. The ratio of the proton to neutron magnetic form factors G^{Q'^)/G^{Q'^) ver- 
sus Q^. The shaded area corresponds to a ±15% deviation from the dipole-fit expectation 
G\f{Q'^)/G^{Q'^) ~ fJ-p/Hn — —1.46 (long-dashed hne). The dot-dashed, dotted, dashed and 
sohd hnes are as in Fig. 1. 
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